Non-uniform chiral phase studied within the Polyakov NJL model 
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We consider how does the introduction of a Polyakov loop affects the spatially inhomogeneous 
quark condensate. The primary result of our work is that the existence of the spatially non-uniform 
chiral phase is confirmed within the Polyakov NJL model in a chiral limit. These findings are 
obtained both in a 3d-cutoff and in a Schwinger (proper time) regularization schemes. 

I. INTRODUCTION 

("^ , The QCD phase diagram is one of the most interesting but still insufficiently explored topics in particle physics. 
p\| ' Both theoretical and experimental efforts are currently taken to get a better understanding of the phase diagram of 
the strongly interacting matter. One of them are studies on the the spatially non-uniform chiral condensate. Since 
j^ ' the existence of the spatially non- uniform chiral phase in a dense quark matter was proposed [l| , many studies have 
been devoted to analyzing this subject [2, 3]. People were interested in determining the range of baryon densities and 
temperatures for which the non-uniform chiral condensate is a ground state of a system. The Nambu - Jona-Lasinio 
model [J], as an effective theory, that describes chiral symmetry breaking was widely used to investigate the QCD 
C^ ■ phase diagram, and the results of the NJL model have been tested due to the change in regularization scheme and 
parameters choice. Only recently, the non-zero quark mass was included and the spataially non-uniform phase was 
studied not only in a chiral limit 0-0] • Nonetheless NJL model appeared to be very useful in analyzing a QCD phase 
rS i' diagram, it has a clear deficiency. NJL theory does not describes the significant aspect of QCD, what is a quark 
confinement. Hopefully, the Polyakov-NJL model was proposed as an extension of the original NJL theory 8, 9(|. In 
t^i the Polyakov-NJL model the background gluon field is introduced and the expectation value of the Polyakov loop can 
r^ ' be an indicator of a transition from a confined into a deconfined regions [10|. In the present paper, we study in the /x 
I— '. - T plane, how does the introduction of a Polyakov loop affects the spatially inhomogeneous condensate. 






II. MODEL 
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^^D The initial lagrangian density of the PNJL model with two massless quarks has a form 

g : H= j ^{i^^D, + /170 ~ m)^j + G [{Ij^^jf + [i^iT^r^f] , (1) 
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where V' is the quark field and /i is the quark chemical potential. The color, flavor and spinor indices are suppressed. 
The vector t is the isospin vector of Pauli matrices. The integration / = /q dT j (fix, where /3 is the inverse 
temperature and the covariant derivative operator D^ = d^ — iA^{x), where A^{x) is a SU(3) gauge field. The 
coupling constant G describes interaction that is responsible for the creation of a quark-antiquark condensate. There 
is also an additional parameter A which defines the energy scale below which the effective theory applies. It is 
introduced through the regularization procedure. We are working in the mean field approximation within the ansatz 
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(^^) = -—cosg-f, {i}}i-^zT''ip) ^ ~—5a3S\TLq-x (2) 

which describes three possible phases: the chiral uniform phase Ch {q — Q,M ^ 0), the non- uniform chiral phase 
NCh (g 7^ 0, M 7^ 0) and the quark matter phase QM {q = 0, M = 0). Because the PNJL model is described in detail 
in a vast literature on this subject [12Hl4| . we only notice some basic assumptions. In our paper we use an effective 
potential of the form [iJl 
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62(T) =6.75- 1.95 ('^]+ 2.625 ('^] - 7.44 ('^') , 63 = 0.75, 64 = 7.5, (3) 



where $ and $ are the traced Polyakov loop and its conjugate. Following Ref. '[15] we examine two scenarios for the 
parameter Tq. First, when Tq is constatnt and equal 208 MeV. Second, when the additional dependence of a Polyakov 
critical temperature Tq on a finite chemical potential is introduced. 
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Next, we assume that the spatial term of a gluon gauge field is vanishing and that the temporal component Aq is 
constant. With above assumptions, the traced Polyakov loop $ is expressed by the gluon filed Aq in the form 



$ = rrL/iV„ L^exp{-l3Ao) 



(5) 



Because we stay in a mean field limit, the expectation values < $ > and < $ > are equal to each other and real [l| 
To determine the ground state of the system we calculate the thermodynamic potential 
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Using the standard method one can obtain an explicite form of the potential ilg 
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The PNJL model is a nonrenormalizable theory, that is why, we apply a regularization procedure that finally defines 
the model. We use a 3d cutoff scheme and a Schwinger (proper time) scheme. Parameters of the model are uniquely 
fixed to reproduce the values of the physical quantities in the vacuum, the pion decay constant /^ = 93 MeV and the 
quark condensate density < uu >=< dd >= —{250MeV)^ [13, llH. After these procedure, thermodynamic potential 
r^o is a fuction of a chiral order parameter M, an absolute value of the wave vector q, a traced Polyakov loop <&, a 
quark chemical potential /z and temperature T. For constant fi and T we numerically solve the mean field equations 
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III. RESULTS 



A. Phase Diagram 



Our aim is to resolve whether in the framework of the PLNJ model, the spatially non-uniform condensate could 
exists and indicate in the fJ. - T plane the region where possibly the NCh phase is a global minimum. We also show, 
how the particular parameters of the thermodynamic potential f^o depend on fi and T. We begin this analysis with a 
brief introduction to the appearance of the phase diagram of the NJL model (with two massless quarks) [l^. In Fig. 
1 we can see the illustration of a ^ - T phase diagram with pre- fixed q = (left panel) and with q j^ (right panel). 
Note that introducing a possibility of an existence of the spatially inhomogencous phase results in a distinct change 
of a critical chemical potential (transition between chiral and QM phase). The point of a chiral phase transition (at 
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FIG. 1: The phase diagrams of a NJL model (with two massless quarks) in a 3d cut-off regularization scheme, GA'^=2.14, 
A =653 MeV. The critical end point (black dot) on a left panel is located at /i =0.28 GeV, T =0.08 GeV. 
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FIG. 2: The phase diagrams of a PNJL model (with two massless quarks) in a 3d cut-off regularization scheme, GA'^=2.14, 
A =653 MeV, To =208 MeV. The critical end point (black dot) on a left panel is located at /i =0.27 GeV, T =0.14 GeV. 



zero temperature) is shifted about 125 MeV towards the larger baryon density. The line of a transition from the Ch 
into the NCh region (right panel) is slightly moved to the left in comparison with a transition line from the Ch into 
the QM phase (left panel). The phase transitions between Ch/NCh, NCh/QM are first order but the transition line 
between Ch/QM (right panel of Fig. 1) is a smooth crossover. The phase transition between Ch/QM (left panel of 
Fig. 1) is of the first order untill the critical point (marked by the black dot). 

Now, we consider the same diagrams but obtained in the Polyakov NJL model. In Fig. 2, there are presented phase 
diagrams, without NCh condensate (left panel) and with NCh phase (right panel). As one can expected, because 
the Polyakov potential is proportional to the fourth power of temperature, the introduction of an effective potential 
for $, affects diagrams mainly in a large temperature region, while at T = 0, the critical chemical potentials remain 
unchanged. As it turns, there is a region where the ground state of the system is realized by the spatially non-uniform 
condensate (gray area). The nature of phase transitions in Fig. 2 is the same as previously. On the left panel of Fig. 
2, we see that temperature of a chiral crossover at /i = is now Tc = 221 MeV. It is a consequence of an interplay 
between the Polyakov effective potential W($, $, T) and the rest part of the thermodynamic potential Jig proportional 
to temperature. The chiral crossover and the deconfinement transition are distinct phenomena, however they are 
linked with each other. In the presence of dynamical quarks, the expectation value of a traced Polyakov loop plays 
a role of an indicator of a deconfinement (not an exact order parameter). We can indicate a region of a very rapid 
change in $ and linked that region with a deconfinement transition. In a scenario when the parameter Tg of a Polyakov 
potential U is constant and equal 208 MeV, the coincidence between a chiral and a deconfinement transitions is valid 
only for small chemical potentials and large temperatures. Apart from these region, the deconfinement transition 



line goes after the chiral one, what is an unphyssical efect 15[ and that is why we do not present the deconfinement 
line in Fig. 2. The partial solution to that situation is suggested in Ref ^^. Additional /i-dependence of a critical 



temperature Tq is introduced ((4]). The phase diagram obtained in this latter scenario is presented in Fig. 3. To(/i) 
is getting smaher with growing chemical potentaial, what results in a fact that a deconfined transition (marked by 
the dashed line) precedes a chiral transition and this relation is valid untill the region of very large fi and small 
temperature. 
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FIG. 3: The phase diagrams of a PNJL model (with two massless quarks) in a 3d cut-off regularization scheme, GA^=2.14, 
A =653 MeV, To = To{fi). The critical end point (black dot) on a left panel is located at fi =0.288 GeV, T =0.085 GeV. Dotted 
line denotes the deconfinement transition. 

The outcomes of the PNJL model are strongly dependent on the choice of parameters and regularization procedure. 
Therefore, we also examine the results obtained in a Schwinger regularization. In Fig. 4, the phase diagram obtained 
in a proper time scheme is presented. To =208 MeV on the left panel of Fig. 4 and on the right, Tq is a function of 
fi. The values of a critical chemical potential and temperature are different than in a 3d cutoff scheme. The chiral 
crossover and the deconfinement transition at zero baryon density almost coincide, while in a 3d cutoff scheme, the 
distance between them was about 20 MeV. Despite these quantitative differences, the main result is common and the 
region of an existence of the NCh condensate is marked in gray. 
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FIG. 4: The phase diagrams of a PNJL model (with two massless quarks) in a Schwinger regularization scheme, GA^=3.78, 
A =1086 IVIeV. On the left panel To =208 MeV, on the right panel To = To{fi). Dotted line denotes the deconfinement 
transition. 



B. Parameter Dependences 



A nice way to show the behavior of the various parameters of the thermodynamic potential VLq is to perform several 
cross-sections through the full phase diagram. In Fig. 5, we present the dependence of a quark constituent mass 
M, an absolute value of a wave vector q and a traced Polyakov loop $ on fi, for respectively T = 20 MeV (left 
panel) and T = 120 MeV (right panel). We see, that this dependence is very similar to that obtained in the original 
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FIG. 5: The dependence of M (gray line), q (open boxes) and "1> (dotted line) on [i, at one plot in a 3d cut-ofT regularization 
scheme, GA^=2.14, A =653 MeV, To =208 MeV. Left panel T = 20 MeV, right panel T = 120 MeV. 
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FIG. 6: The dependence of M (gray line), q (open boxes) and $ (dotted line) on T at one plot in a 3d cut-off regularization 
scheme, GA^=2.14, A =653 MeV, To =208 MeV. Left panel /^ = MeV, right panel p = 260 MeV. 



NJL model [1^, [T^ . Both a transition from the uniform chiral to the non- uniform chiral and a transition from the 
non-uniform phase to the quark matter phase are of the first order. The value of $ is gradually growing with an 
increase of a chemical potential. It is difficult to find any physical interpretation for a discontinuity of the $ value 
in a chemical potentials corresponding to Ch/NCh and NCh/QM transitions. Nevertheless, we remember that in 
our case, the traced Polyakov loop $ is only an indicator of a deconfinement, and so we are interested mostly in the 
qualitative results. In Fig. 6, the additional analysis is presented. The dependence of M, q and $ on temperature 
is examined for a constant chemical potential. On the left panel of Fig. 6, one can observe a mentioned earlier 
merging of a chiral crossover and a deconfinement transition. Region of a very rapid change in $ coincides with the 
disappearance of M. On the right panel of Fig. 6 (/i = 260 MeV) the growth of $ is more gradual. Again, there are 
discontinuities of a $ value, one associated to the Ch/NCh transition and the other to the NCh/QM transition. On 
the left panels of Fig. 5, 6, (scenario when To = 208 MeV) it is visible that while the coincidence between a chiral 
and a deconfinement transitions makes sense for small chemical potentials and high temperatures, in the rest part of 
a diagram, this agreement is not entitled. Introduction of a /^-dependence of a Polyakov critical temperature causes 
that the region of a rapid growth in $ appears for smaller chemical potentials and the growth is much steeper. It is 
visible in Fig. 7, where the M, q and $ values are plotted for respectively T = 120 MeV left panel and fi =260 MeV 
right panel. 
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FIG. 7: The dependence of M (gray line), q (open boxes) and $ (dotted line) at one plot in a 3d cut-off regularization scheme, 
GA2=2.14, a =653 MeV, To = To(At). Left panel T = 120 MeV, right panel p = 260 MeV. 

IV. CONCLUSIONS 

The primary result of our work is that the existence of the spatially non-uniform chiral phase is confirmed within 
the Polyakov NJL model in a chiral limit. These findings are obtained both in a 3d cutoff regularization and 
in a Schwinger (proper time) schemes. We have examined two scenarios of the critical Polyakov temperature Tq 
dependence. One, when Tq is constant and equal 208 MeV, the other, when the Tg value is lowering with growing 
/z. For both of these scenarios, in comparison to the original NJL model, the temperature of a chiral crossover in 
/J, = changes from 190 to 221 MeV (3d cutoff) and from 152 to 189 MeV (proper time), while the critical chemical 
potentials at T = between Ch/NCh and NCh/QM remain unchanged. As a consequence the line of a transition 
from chiral phases into the quark matter phase is always above the same line obtained in the original NJL model. 
This effect is stronger for constant Tq. Introduction of a Tq{^) dependence causes a better agreement of a chiral and 
a deconfinement transitions. The future work that should be done to get a better understanding of the PNJL phase 
diagram is to concider simultaneously the non-zero current quark mass and inhomogeneous condensate. In the present 
paper we stay in a mean field approximation, the analysis of a diagram away from this limit is also worth consideration. 
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